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Abstract 

We propose, for the computation of the Higgs mass spectrum and cou- 
phngs, a renormahzation-group improved leading-log approximation, where 
the renormalization scale is fixed to the top-quark pole mass. For the case 
niA ~ MsusY, our leading- log approximation differs by less than 2 GeV 
from previous results on the Higgs mass computed using a nearly scale in- 
dependent renormalization-group improved effective potential up to next- 
to- leading order. Moreover, for the general case rriA ~ ^susy? we provide 
analytical formulae (including two-loop leading-log corrections) for all the 
masses and couplings in the Higgs sector. For Msusy ~ 1-5 TeV and arbi- 
trary values of rriA, tan/? and the stop mixing parameters, they reproduce 
the numerical renormalization-group improved leading-log result for the 
Higgs masses with an error of less than 3 GeV. For the Higgs couplings, our 
analytical formulae reproduce the numerical results equally well. Compar- 
ison with other methods is also performed. 
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1 Introduction 



The most appealing extension of the Standard Model (SM) is the Minimal Supersym- 
metric Standard Model (MSSM), which can provide a technical explanation for the 
hierarchy stability from Mpi to the electroweak scale. From the experimental point 
of view, the MSSM is also attractive and presents clear and distinct signatures. In 
particular, its Higgs sector predicts the existence of a light CP-even state with a mass 
which cannot exceed a value of ~ 150 GeV, unlike the case of the SM Higgs whose mass 
can reach much larger values (~ 1 TeV). The precise bounds on the mass spectrum of 
the MSSM Higgs bosons strongly depend on the value of the top-quark mass Mt lH, @]. 

In view of Higgs searches at future colliders, and of the fact that the top-quark 
mass is being measured at the Tevatron [|, it is of the highest interest to provide 
determinations as precise as possible of the radiatively corrected Higgs mass spectrum 
and Higgs couplings and, in particular, to find accurate upper bounds on the mass of 
the lightest Higgs boson. To this end one can assume that Mw <^ Msusy ~ O(few) 
TeV, where Msusy is the characteristic supersymmetry particle mass scale and its 
upper limit is determined by naturalness arguments. In this case the upper bound 
on the lightest Higgs mass is reached when the mass of the CP-odd Higgs, m^, is of 
the order of Msusy, and the theory below Msusy is the SM with threshold effects at 
-^susY- For values of < Msusy the theory below Msusy is the two-Higgs doublet 
model, and then, two CP-even and one CP-odd neutral Higgs bosons, and two charged 
Higgs bosons appear in the physical spectrum below Msusy- 

It is the aim of this letter to explore the whole range < Msusy and provide an 
accurate numerical evaluation as well as analytical formulae for the whole Higgs mass 
spectrum and Higgs couplings of the MSSM. We shall also compute the range of vahdity 
of our approximations, as well as make comparison with other methods currently used 
in the literature. 

2 The case ttia ~ M^gusY 

This case is well described, as we mentioned above, by the SM with the quartic coupling 



In the above, g2 and gi are the SU (2) x f/(l)y gauge couplings, tan /3 = f 2/f i is the ratio 
of the vacuum expectation values of the neutral components of the two supersymmetric 
Higgs bosons, is the SM top-quark Yukawa coupling and Xt is the stop mixing 
parameter ^: 
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at Msusy, and one-loop threshold contributions at that scale 
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^In this section we are neglecting the sbottom mixing parameter. 
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At = At — fi cot p. 
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The supersymmetric scale Msusy is taken as the squared root of the arithmetical 
average of the diagonal entries in the squared squark mass matrix and, for simplicity, 
we assume that all supersymmetric particle masses are of order Msusy- A more detailed 
discussion about the definition of Msusy and the possible effect of light sparticles will 
be given in sections 4 and 5. 

The Higgs boson mass was first determined by renormalization-group resummation 
of all- loop leading- log (LL) corrections |, |]. Some next-to- leading- log (NTLL) 
corrections were further introduced |T0| and finally a complete NTLL analysis was 
performed in Refs. ||Tl], |12|. 

In Ref. [|12| physical pole masses for the top-quark and Higgs-boson were used and 
the dependence with respect to the renormalization-group scale in the NTLL approx- 
imation was explicitly analysed. One of the main issues in |T2[ was the comparison 
between the LL and NTLL approximations. As expected, the LL approximation shows 
strong scale dependence, while the NTLL is almost scale independent. This implies 
not only that the choice of the renormalization-group scale is irrelevant when working 
in the NTLL approximation, but also that a judicious choice of the renormalization 
scale in the LL approximation can yield very accurate results. It turns out that, in- 
dependently of the value of the stop mixing parameter, the scale where both results 
coincide is very close to the pole top-quark mass Mt (see e.g. Fig. 8 of Ref. [T^). 

This result can also be understood from Ref. 0] where the physical Higgs mass is 
related to the quartic coupling in the MS scheme, A(t). One can write 
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where 
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is the running Higgs mass, and v = 174.1 GeV. From |jT3|, |T4|, and neglecting scalar 
and electroweak gauge couplings, we can write 
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where Q is the renormalization scale and the function Z{x) is given, for z > 1/4, by 

Z{x) = 2V4a; - 1 tan"^(l/V4a; - 1) . (7) 

We can see from (^) that choosing Q = Mt, the difference between the pole mass and 
the running mass is determined hj e {e < 1.7% for > 60 GeV and Mt < 200 GeV). 
For instance for Mh = 90 GeV and Mt = 180 GeV, then e = 1.3%. Other values of Q 
(e.g. Q = Mz as can be usually found in the literature |T^) can provide much larger 
corrections, depending on the values of M^ and Mt, unless an appropriate decoupling 
of the top-quark effects below the top-quark mass scale is performed. 
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For the Higgs mass, we shall therefore adopt the value provided by (^) with Q = 
and work in the LL approximation as suggested by (H). We compare, in Fig. 1, the 
all-loop RG improved NTLL Q (solid) and LL (dashed) lines for Msusy = 1 TeV. 
We consider the case of zero mixing, Xt = 0, and the case Xt = 6, which maximizes 
the lightest CP-even Higgs mass, large values of tan/5 (tan/5 = 15) and small values 
of tan/3 [the infrared (IR) fixed point solution with sin/? ~ (200 GeV)/Mt]. The RG 
improved LL aproximation involves the one-loop RG running of the couplings, the one- 
loop threshold effects AA, Eq. (H), from the decoupling of the supersymmetric particles 
at MgusY, and the computation of the Higgs masses at the scale Q = Mt. Moreover, 
in both approximations the pole top-quark mass Mt was computed from the on-shell 
running mass mt through the corresponding one-loop QCD correction factor 

Mt 

= l + i«3(M.)- 

We see that, for all values of the top-quark mass, tan P and the stop mixing parameter, 
the solid and dashed curves agree with each other with an accuracy better than 2 GeV. 

We have also worked out an analytical approximation to the numerical all-loop 
renormalization-group improved LL result, including two-loop leading- log effects. It is 
given by 
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where the angle /3 is defined here at the scale Msusy 
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<^log%^. (10) 

and 

""^'"'^ = l + |ia3(M.)log(A/,7A/|)' 

where 63 = 11— 2A^//3 is the one-loop QCD beta function and Nf is the number of quark 
flavours, which is equal to 5 at scales below Mt. Notice that Eq. (D includes the leading 
D-term correction (9(M|m^). As was pointed out in Ref. |1^, this contribution cannot 



be neglected, since, for the experimental Mt range, it can account for a (negative) shift 
as large as ~ 5 GeV. Formula (^ is also plotted in Fig. 1 (dotted lines), for the same 
values of the supersymmetric parameters as before. We can see that the analytical 
approximation reproduces the numerical results within an error of less than 2 GeV in 
all cases. 

Finally, in Fig. 2 we have plotted the Higgs mass, in the three different approxi- 
mations, for Mt = 175 GeV, as a function of Msusy- As expected the NTLL (solid) 
and LL (dashed) curves coincide within less than 2 GeV difference. The analytical 
approximation also remains within this level of accuracy if Msusy ~ 1-5 TeV. Hence, 
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it follows that for higher values of Msusy, agreement with numerical results would 
demand the addition of 0{t^) terms (three-loop leading- log) , but for the theoretically 
preferred values of the supersymmetry breaking scale, Msusy ~ 1 TeV, the analytical 
formula provides excellent results. For values of Msusy < 500 GeV, the maximal value 
of the mixing parameter Xt is restricted (see discussion in section 4). 



3 The case rriA ~ AfgusY 

In the more general case, < Msusy, the effective theory below the scale Msusy H] is 
a two-Higgs doublet model where the tree level quartic couplings can be written in terms 
of dimensionless parameters Aj, z = 1, . . . , 7, whose tree level values are functions of the 
gauge couplings, and with one- loop threshold corrections AA,, i = 1, . . . , 7, expressed 
as functions of the supersymmetric Higgs mass fi and the soft supersymmetry breaking 
parameters At, and rriA- We are using the conventions of Ref. Motivated by 
the previous agreement among the NTLL and LL approximations and the analytical 
formulae, we can safely extrapolate the analytical approach, to obtain the two-loop 
results for the Higgs spectrum within the MSSM. 

The CP-even light and heavy Higgs masses and the charged Higgs mass are given as 
functions of tan/3, Msusy, At, Ab, fi, the CP-odd Higgs mass rriA and the physical top- 
quark mass Mt related to the on-shell running mass rrit through Eq. (^. We present 
below the masses and mixing angle of the Higgs sector as functions of the parameters 
Aj, for which we have derived the corresponding analytical formulae and which are 
the equivalent in the two-Higgs doublet model of the analytical approximation (^) for 
the case of the SM. For completeness, we include the bottom Yukawa coupling effects, 
which may become large for values of tan/5 ~ mt/mi,, where is the running bottom 
mass at the scale Mt. 

The two CP-even and the charged Higgs masses read 



where 
with 



ml^H) = ' 2 ' ^^^^ 



TrM2 ^ y^(TrM2)2 -4det 
2 

m^± = m\ + (As - A4)^;^ (13) 
TrM^ = Ml^ + MI2 ; det M^ = M^^M^ - (mI^'^ , (14) 



M^2 = 2f ^ [sin (3 cos /5( A3 + A4) + Ae cos^ + A7 sin^ /?] — sin (3 cos (3 
Ml^ = 2t;2[Aicos^/3 + 2A6Cos/3sin/5 + Assin^/?] +m^sin2/3 (15) 
= [A2 sin^ /5 + 2 A7 cos /3 sin (3 + X5 cos^ (3] + m\ cos^ (3. 

The mixing angle a is equally determined by 

9 A/f2 

sin 2a = (16) 

'(TrM2)2 -4detM2 
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The above quartic couplings are given by 
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which contain the same kind of corrections as Eq. 
contributions, and where we have defined, 



including the leading D-term 
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All quantities in the approximate formulae are defined at the scale Mt. In particular 
a^lMt) is given in Eq. ([Til) , and 



ht = mt{Mt)/{v sin p) 
hh = mh{Mt) / {v cos P) 



(26) 



are the top and bottom Yukawa couplings in the two-Higgs doublet model. 

For rriA < Mt, tan/5 is fixed at the scale m^, while for 171^ > Mt, tan/5 is given 

by [g 



tan P{Mt) = tan/?(myi) 



(27) 



For the case in which the CP-odd Higgs mass tua is lower than Msusy, but still 
larger than the top-quark mass scale, we decouple, in the numerical computations, the 
heavy Higgs doublet and define an effective quartic coupling for the light Higgs, which 
is related to the running Higgs mass at the scale through 



A(mA) 



mh[mA) 
2v^ 
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The low energy value of the quartic coupling is then obtained by running the SM 
renormalization-group equations from the scale down to the scale Mt. In the 
analytical approximation, for simplicity we ignore the effect of decoupling the heavy 
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Higgs doublet at an intermediate scale. We partially compensate the effect of this 
approximation by relating the value of tan/5 at the scale Mt with its corresponding 
value at the scale niA through its renormalization-group running, Eq. (pTf). 

In Figs. 3, 4 and 5 we plot mh, mn and m^, respectively, as functions of m^, 
for Mt = 175 GeV and the four cases corresponding to tan/3 = 15 and tan/? ~ 1.6 
(which is the lowest possible value of tan /? for this Mt and corresponds to the IR 
fixed point solution) and for large and zero mixing. Solid curves correspond to the 
numerical renormalization-group improved LL approximation and dashed curves to 
the analytical approximation depicted in the previous formulae. For the charged Higgs 
mass, we have also plotted the case fi = 2MgusY and tan/? = 15, since, unless large 
values of tan/3 ~ rrit/mi, are considered, its radiative corrections depend on /x but not 
on At or A), (see Eqs. ( pT] ) and (p2D). We can observe that, in all cases, the agreement 
between the Higgs masses in the numerical renormalization-group improved calculation 
and the analytical approximation is, either of the order of, or better than, 3 GeV. 



Vertex 


Coupling 


{h,H}W,W, 
{h,H}Z^Z, 

{h, H, A}uu 

{h, H, A}dd 

{h,H}AZ^ 


ig2Mwg^iu{sin{P - a), cos(a + P)} 
192 5'Mi'{sin(/5 a),cos(a + /3)} 

COS^ dw 

oU^ • ^|cosa,sma, i-f^cos/S} 

o «i sma,cosa, «75sm/?} 
2 \Mw/ cosp 

'^^/^^^ {cos(/? -a),- sin(/5 - a)} 
2 cos Ow sm Ow 



Notice that the knowledge of the radiatively corrected quart ic couplings Xi, i = 
1, . . . , 7, and hence of the corresponding value of the Higgs mixing angle a, permits the 
evaluation of all radiatively corrected Higgs couplings. For instance some important 
Yukawa and gauge Higgs couplings Ijl^ are listed in the table above, where {k^) 
is the incoming (outcoming) CP-odd (CP-even) Higgs momentum. Furthermore, the 
trilinear Higgs couplings can be explicitly written as functions of Aj, i = 1, . . . , 7, a 
and /3 P, |18]. To check the goodness of our analytical approximation on the different 
couplings we plot in Fig. 6, for the same four cases as in Figs. 3-5, cos a as a function 
of tua, for Mt = 175 GeV. As can be seen from Fig. 6, we find excellent agreement 
between the values obtained with the numerical and analytical approximations. 
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4 On the expansion variables 



In this section we shall comment on the definition of MgusY- As we mentioned above, 
the scale Msusy, appearing in the quartic coupling expressions, should be associated 
with a characteristic squark mass scale. The renormalization group analysis considered 
in this work is based on the assumption of a single step decoupling of the squark fields, 
which is valid only if the mass splitting among the squark mass eigenstates is small. 
More quantitatively, the validity of our expansion requires, 

2 2 

— mf 

^ < 0.5, (29) 

mi + m~ 

where m?^ and m?_^ are the stop squared mass eigenvalues. This restriction also applies 
to all existing Higgs mass analyses at the next-to leading order [10-12]. Observe that 
Eq. (p9D automatically ensures the absence of tachyons in the stop sector. Furthermore, 
under the validity of Eq. (|29|), the scale Mgugy may be safely defined as the arithmetic 
average of the stop squared mass eigenvalues, 

MlusY = '\ (30) 

Observe that our work, as well as all other RG Higgs mass analyses performed in the 
literature, also relies on an expansion of the effective Higgs potential up to operators 
of dimension four. The contribution of higher dimensional operators may be safely 
neglected only if 2\MtAt\ < M^^^y and 2\Mtfx\ < M|usy- 

Strictly speaking, the operator expansion was performed in the symmetric phase, 
where the squared of the energy scale at which the stops decouple is not the average 
of the stop squared masses, Msusy, but the average of the stop soft supersymmetry 
breaking squared mass parameters, = (mQ + m^)/2. Higher dimensional operators 
in the symmetric phase can only be neglected if M|usy ^ ^t- ^^^^ range of 
parameters, the distinction between Msusy and the real decoupling scale has no impact 
on our approximations. In general, our expressions give a reliable approximation to the 
Higgs masses whenever the above restriction on the stop mixing parameter, Eq. (pQ]), 
holds, and Msusy ~ SMj. Moreover, for the large tan/? case, for which the bottom- 
Yukawa contribution becomes relevant, we have assumed that the sbottom masses are 
of order of the stop ones and that similar bounds on their mixing mass parameters are 
fulfilled. 

Physically it is clear that, in the region where higher order operators are unsup- 
pressed, Msusy ~ 2Mi, the leading logarithm expansion of the Higgs masses should 
contain powers of log(MgusY/^t^) rather than log(m^/Mj^). This should be the case, 
since, in the supersymmetric limit for the top/stop sector, Msusy = Mt, and zero 
mixing, one should recover the tree level values for the Higgs masses. One may ask 
how accurately our analytical expressions may determine the Higgs masses for values 
of Msusy ~ 2Mt, for which, strictly speaking, one is beyond the limits of the present 
expansion. In fact, in order to check the reliability of our approximation for the case 
in which light stops, with masses of the order of the top mass, are present, an al- 
ternative computation, valid in that regime, should be used. The expressions for the 



8 



CP-even Higgs mass matrix elements computed in Ref . , and based on the one-loop 
effective potential in the MSSM (keeping the leading supersymmetric corrections from 
the stop/sbottom sectors) can be useful for the comparison. Observe that the MSSM 
effective potential is valid only at scales of the order of, or larger than, the charac- 
teristic stop supersymmetry breaking terms. One cannot consider the MSSM effective 
potential at scales below the squark masses since, at those scales, the squarks are al- 
ready decoupled, possibly leaving threshold corrections to match the effective potential 
beyond and below the supersymmetric threshold. 

Once recognized that the MSSM effective potential is valid at the scale defined 
by the squark masses, one can try to improve its performance by means of its RG 
improvement [|l^. Were the RG improved effective potential exactly scale independent, 
one could minimize it at any scale (e.g. at Msusy), defining all the running parameters 
on which the effective potential depends also at that scale, and rewrite them, in turn, 
as functions of the corresponding parameters defined at the scale Mt. Then, taking 
advantage of the scale independence of the effective potential one could run the Higgs 
masses from the high scale to Mt using the anomalous dimensions of the Higgs fields. 
However, the RG improved one-loop effective potential is not exactly scale independent 
(even after including all LL and NTLL corrections) and hence, when minimizing it at 
MsusY (see e.g. Figs. 1 and 3 in Ref. [jl2|) one is implicitly including an error in 
the lightest CP-even Higgs mass determination. This error is very small for values of 
-^susY of the order of the weak scale, but becomes significant for larger scales (e.g. 
MsusY ^ 500 GeV). 

We therefore consider the RG improved MSSM effective potential, which is min- 
imized at the scale Msusy, at which we determine the neutral Higgs boson squared 
mass matrix. All the parameters in the mass matrix are running parameters, defined 
at the scale Msusy (including the vacuum expectation values Vi and V2). They are 
related to parameters defined at the scale Mf by the corresponding evolution with the 
/3- and 7-functions, that we compute to one-loop LL (i.e. first order in the expansion 
parameter log(MgusY/^t^)) and keeping only the Yukawa and strong gauge couplings, 
as in our analytical expressions of sections 2 and 3. We then evolve the squared mass 
matrix, with the corresponding anomalous dimension matrix, from MgusY to M^. In 
Fig. 7 we plot the lightest CP-even Higgs mass obtained from our analytical expressions 
(solid lines) and that obtained from the RG improved effective potential of the MSSM 
(dashed lines), for a pole top-quark mass Mt = 175 GeV, vanishing mixing At = fi = 0, 
iTiA = Msusy (to simplify the comparison), and large (tan/3 = 15) and small (IR) 
values of tan We can see that the agreement between the Higgs masses obtained by 
using our analytical expressions and the RG improved MSSM effective potential is very 
good for all values of Msusy ~ 500 GeV, and in particular for Msusy ~ 2Mt, for which 
the latter should give a reliable description. Hence, we see that, independently of the 
value of tan /?, our analytical expressions give a very good approximation to the Higgs 
mass values, in the case of no left-right stop mixing, even in the regime Msusy ~ 2Mj. 
For the case of stop mixing considered in Fig. 8, observe that, independently of the 
value of tan /?, the results are also in very good agreement for values of Msusy ~ 500 
GeV. Similar conclusions may be obtained for any value of the stop mixing parameter. 

Observe that, in order to compute pole Higgs masses, Higgs vacuum polarization 
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contributions should be added to our computation. As we discussed in section 2, 
by fixing the renormahzation scale to the pole top-quark mass, we make the top- 
dependent vacuum polarization contributions small. The stop-dependent contributions 
are renormalization-scale independent. For light stops, these contributions are opposite 
in sign to the top leading ones and of the same order of magnitude; they become very 
small if the stops are heavy. In general, for our choice of the renormahzation scale, 
the vacuum polarization contributions are small, and would only be required in a 
computation of Higgs masses going beyond the approximation presented in the present 
work. 

Summarizing, even for stop masses of order of the top mass, our analytical expres- 
sions give a very good approximation to the Higgs mass value. Observe that, since 
we have applied one step decoupling on the stops, the constraint on the stop mass 
splitting, Eq. (^9]), should always be held. Finally, in Figs. 7 and 8 we also plot, for 
comparison, the values obtained from the MSSM effective potential while ignoring the 
RG improvement, and evaluating all running masses and couplings directly at the scale 
M-t (dotted lines). In particular, we use, for the top mass, the running value mt{Mt). 
This approximation overestimates the Higgs mass for basically all values of Msusy- 



5 Conclusions 

In this work we present a one-loop improved renormalization-group analysis of the 
lightest CP-even Higgs mass, which, for the case rriA ~ AfgusY, reproduces the existing 
two-loop renormalization-group improved effective potential results within an error of 
less than 2 GeV. For the general case < Mgusy, we provide an analytical approxi- 
mation (including two-loop leading-log corrections) for the effective quartic couplings, 
which determines the Higgs mass spectrum within an error of less than 3 GeV with 
respect to the values obtained using the improved one loop RG approximation. Our for- 
mulae are also useful to write radiatively corrected expressions for the Higgs couplings, 
which are of interest in computing physical processes. 

Our analysis takes into account the dominant top and bottom quark Yukawa (as well 
as strong gauge) coupling effects, by assuming a characteristic supersymmetric scale, 
-^susY; and including squark mixing effects. In some realistic scenarios the chargino 
and neutralino masses may become much smaller than the characteristic squark mass 
scale. Light charginos and neutralinos have a small and negative effect on the lightest 
CP-even Higgs mass, of order 2-3 GeV [0, which may be added to our analytical 
approximation at the one-loop order. 

Our approximation relies on a one step decoupling of all supersymmetric particles 
and, in particular, on the stop mixing parameter requirements discussed in section 4. 
A computation of the Higgs masses and couplings going beyond the approximation 
presented in this work, would demand the derivation of the full two-loop RG improve- 
ment of the one-loop effective potential, and polarization corrections to running Higgs 
masses, including the squark and heavy Higgs contributions, the analogous to what 



has been done for the case of large squark and CP-odd Higgs masses in Ref. |jT2|. We 
expect to come back to this subject elsewhere. 
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Figure 1: The lightest CP-even Higgs mass as a function of the physical top-quark mass, 
for MsusY = 1 TeV, evaluated in the limit of large CP-odd Higgs mass, as obtained 
from the two- loop renormalization-group improved effective potential (solid lines), the 
one-loop improved RG evolution (dashed lines) and the analytical formulae, Eq. (^ 
(dotted lines). The four sets of lines correspond to: a) tan f3 = 15 with maximal squark 
mixing, Xt = 6; b) tan /? = 15 with zero mixing, Xt = 0; c) the minimal value of tan f3 
allowed by perturbativity constraints for the given value of Mj (IR fixed point), with 
Xf = 6; and, d) tan (3 the same as in (c) with zero mixing. 
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Figure 2: The lightest CP-even Higgs mass as a function of the super symmetric scale 
-^susY, for Mt = 175 GeV, and evaluated in the limit of large CP-odd Higgs mass. 
Solid, dashed and dotted lines are as in Fig. 1. The different sets of curves correspond 
to the values of Xt and tan/3 for cases (a) to (d) of Fig. 1. 
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Figure 3: The lightest CP-even Higgs mass as a function of the CP-odd Higgs mass 
for a physical top-quark mass = 175 GeV and Msusy = 1 TeV, as obtained from 
the one-loop improved RG evolution (solid lines) and the analytical formulae, Eq. ([I2D, 
(dashed lines). The four sets of curves correspond to: a) tan/? = 15 with large squark 
mixing, = 6 (/i = 0); b) tan/3 = 15 with zero mixing Xt = ^ = 0; c) the minimal 
value of tan f3 allowed by perturbativity constraints for the given value of Mt (IR fixed 
point), tan/3 ~ 1.6, with large squark mixing; and, d) tan/3 ~ 1.6 with zero mixing. 
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Figure 4: The heaviest CP-even Higgs mass as a function of the CP-odd Higgs mass 
for the same set of parameters as in Fig. 3. 
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Figure 5: The charged Higgs mass as a function of the CP-odd Higgs mass for the same 
set of parameters as in Fig. 3 (upper overlapping curves), and for fi = 2Msusy and 
tan/5 = 15, for which the radiative corrections become observable (solid and dashed 
lower curves). 



17 



100 200 300 400 500 

[GeV] 



Figure 6: Plot of cos a, as a function of the CP-odd Higgs mass for the same set of 
parameters as in Fig. 3. 
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Figure 7: Comparison between the results for the Higgs mass obtained with the analyt- 
ical approximation (solid hnes), the RG improved one- loop MSSM effective potential 
(dashed lines) and the unimproved one-loop MSSM effective potential (dotted lines), 
as described in section 4, for Mt = 175 GeV, itla = Msusy, fnq = and At — ji — Q. 
The lower set corresponds to tan/3 = 1.6, and the upper set to tan/3 = 15. 
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Figure 8: The same as in Fig. 7 but for — Msusy- 
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